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ABSTRACT 

We present images and initial results from our extensive Spitzer Space Tele- 
scope imaging survey of the W5 Hll region with the Infrared Array Camera 
(IRAC) and Multiband Imaging Photometer for Spitzer (MIPS). We detect dense 
clusters of stars, centered on the O stars: HD 18326, BD +60 586, HD 17505 
and HD 17520. At 24 /xm, substantial extended emission is visible, presumably 
from heated dust grains that survive in the strongly ionizing environment of the 
Hll region. With photometry of more than 18000 point sources, we analyze 
the clustering properties of objects classified as young stars by their IR spectral 
energy distributions (a total of 2064 sources) across the region using a minimal- 
spanning-tree algorithm. We find ~40-70% of infrared excess sources belong to 
clusters with >10 members. We find that within the evacuated cavities of the 
H II regions that make up W5, the ratio of Class II to Class I sources is ~7 times 
higher than for objects coincident with molecular gas as traced by ^^CO emission 
and near-IR extinction maps. We attribute this contrast to an age difference 
between the two locations, and postulate that at least two distinct generations 
of star formation are visible across W5. Our preliminary analysis shows that 
triggering is a plausible mechanism to explain the multiple generations of star 
formation in W5, and merits further investigation. 
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type — stars: formation — stars: pre-main sequence 
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Introduction 



Star formation is a self-regulating process — once massive stars form they immediately 
begin to disrupt their natal environment with their stellar winds and the emission of ionizing 
radiation. Eventually their parental molecular clouds are destroyed, halting further star 
formation. However, it has also been argued that the energy input by these massive stars 
can promote and induce subsequent star formation in the surrounding molecular gas before 
it disperses, above that which wou ld be produced without external forcing. This process 



is given the name 'triggering' (see lElmegreenl Il998l : IZinnecker fc Yorkd 120071 . for reviews 



of this subject). It is vital to understand the balance of cloud destruction and triggered 
star formation if we are to develop a theory that explains the morphology and evolution 
of star forming regions, star formation efficiencies and the initial mass function of stars in 
clusters. It also has relevance for star formation on galactic scales in understanding how 
star f ormation progresses with time, and indeed how spiral structure of galaxies evolves with 
time (ISeiden fc Schulmanlll990l : iJungwiert &: Palouslll994j ). Feedback in star formation can 
be clearly observed in bright-rimmed clouds, where an edge-on molecular cloud is externally 
illuminated by nearby young massive stars, creating a cross-section of the photo-evaporation 
process as the ionization fronts they produce advance into the molecular cloud. 

Two triggering mechanisms are of interest in regard to W5. The first is the creation of 
an ionized H ll region bubble by an initial generation of massive stars within their parental 
molecular cloud, and its s ubsequent expansion. In this 'collect and c ollapse' me c hanisra , 
investigated analytically by IWhitworth et al.l (Il994j ) and numerically by iDale et al.l (l2007al ) , 
the expansion creates a shock front that sweeps up neutral material ahead of it as it progresses 
outward. The gas accumulated eventually exceeds a critical threshold for collapse and gives 
rise to a second generation of star formation. Secondly, inhomogeneity in the ISM often leads 
to small clumps of material remaining exposed inside an H ll region. In this environment, the 
high pressure of the surroun ding ionized gas ha s been suggested as a mechanism to compress 
the clumps and form stars (jStutzki et al.lll988l ). 



^The star forming region W5 is a part of the chain of mo 

19581 ). Its distance has not been definitively establis hed — 



ecular clouds W3/4/5 (IWesterhout 



Becker fc FenkartI (Il97ll ) found a 



photometric distance of 2.2 kpc, but more recently, iHillwig et al.l (120061 ) found that a dis- 
tance o f 1.9 kpc gave t he most consistent results of evolutionary model fits to stellar radii 
in W5. IXu et al.l (120061 ) found a distance to the neighboring region W30H of 1.95±0.04 kpc 
using maser parallaxes — we adopt a distance to W5 of 2 kpc as a conservative intermediate 
between these different estimates. W5 is a relatively isolated star forming region, with an 
apparently simple morphology. Optical imaging shows it is made up o f two roughly circular 
adjoining Hll regions W5 East and W5 West (IKarr fc Martini l2003l ). containing one and 
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four O stars respect ively. At least two of the four O stars in W5 West are multiple systems 



( iHillwig et al.ll2006l ). Both ^^CO and 21 cm radio emission (see iNormandeau et al.l Il996l ) 
demonstrate the same overall shape — ^^CO emission traces the molecular hydrogen gas in 
W5, which appears as two broken rings of emission. Figure 1 shows a map of ^^CO (A = 
2.6 mm) emission integrated over the range —28 to —47 km s~^, made using observations 
with the 14m FCRAO telescope. The bulk of the emission is found between a velocity of 
—31 and —47 km s~^, and thus is lik e ly rou ghly in the same plane. In this image, following 
the naming scheme of IWilking et al.l (119841 ). to the northwest of HD 17505 is the molecular 
cloud W5NW, in between HD 18326 and BD +60 586 is W5NE, and to the east of HD 
18326 is W5A. Due to this relatively simple morphology, W5 presents a useful test case for 
investigating models of triggered star formation and the influence of massive star formation 
on its surroundings. 

In this paper we present initial results from our mid-IR Spitzer survey of W5 with 
the IRAC and MIPS instruments. In § 2 we describe our observations and data reduction 
techniques and our classification and clustering analysis in § 3. In § 4 we perform a trial, 
simplified investigation of triggered star formation models as a means to explain the observed 
distribution of young stars in different evolutionary states across the region. In § 5 we present 
our conclusions and directions of future work. 



2. Observations and Data Reduction 



W5 was observed with the Spitzer IRAC instrument (IFazio et al.ll2004l ) in all four bands 
(3.6, 4.5, 5.8 and 8.0 /im). The observations were broken down into three rectangular 
Astronomical Observing Requests (AORs) covering ~1.8xl.6 degrees, in order to observe 
at multiple rotation angles and help minimize artifacts aligned along columns or rows of 
the array. In Table 1 we list the dates and coordinates of each AOR. Each AOR had a 
coverage of 1 High Dynamic Range (HDR) frame. HDR mode results in a 10.4 second 
and 0.4 second exposure being taken at each position in each map. We used software tools 
(clustergrinder) developed by one of us (RAG) to produce final image mosaics from these data 
in each wavelength band. Clustergrinder incorporates all necessary image treatment steps 



for example, saturated pixel processing and distortion corrections (see lGutermuth et al.ll2008 



for a more complete description of the processing performed) . Clustergrinder uses the short 
0.4 sec exposures only in saturated or near-saturated regions, so that the combined map has 
an effective total integration time of 3x10.4 = 31.2 sec in most of the overlapping areas. We 
also incorporated in ou r data processing archival data covering AFGL 4029 (from Spitzer 
GTO program PID 201 1 Allen et al1l2005[ ). at the eastern end of W5. 
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Fig. 1. — ^^CO (1-0) greyscale image of W5 from FCRAO telescope (integrated corrected 
for main beam efficiency rjMB = 0.45, ranging from 0-110 K km s~^), oriented North- up and 



East-l eft. We label molecular clouds in black following the naming scheme of IWilking et al. 
(Il984j ). We label the two Hll regions in white. Prominent O stars are marked with white 
asterisks: HD 18326 in W5 East, and in W5 West from North-East to South- West are BD 
+60 586, HD 237019 and HD 17505 (upper of pair), HD 17520 (lower of pair). 



Table 1. 



AORKEY 


Date 


Coordinates 


IRAC Reduction 




(UT) 


(J2000) 


Pipeline ver. 


14507776 


2006 September 20 


RA 02 53 32.9 


814.0.0 






Decl. +60 25 59.1 




14507008 


2006 September 28 


RA 02 53 43.5 


814.0.0 






Decl. +60 21 57.2 




14508544 


2007 February 16 


RA 02 53 54.2 


815.3.0 






Decl. +60 17 55.2 





- 5 - 



The MIPS ( iRieke et al.ll2004j ) observations were carried out on 2006 February 23 UT 
under our GO-2 program, PID 20300. Images were taken in scan map mode using the medium 
scan speed for an average exposure time of 41.9s per pixel once frames were combined. The 
raw data were processed with pipeline version S13.2.0. We produced final mosaics using 
the MIPS instrument team Data Analysis Tool, which calibrates the data and app lies a 



distortion correction to each individual exposure before combining (IGordon et al.ll2005l ). We 
used only the 24 /im band data for our analysis in this paper, since strong background 
emission dominates at the longer wavelength (70 and 160 /im) bands of MIPS, and lower 
sensitivity reduces the number of detectable objects to an un-useful level for the present 
study. 

In Figure 2 we present a 3-color IRAC image of W5, using the color scheme: blue=3.6 
/im, green=4.5 fim and red=8.0 fim. In Figure 3 we show our MIPS 24 fim mosaic. Figure 
4 presents a composite image incorporating MIPS, with blue=4.5 /im, green=5.8 /im and 
red=24 /xm (MIPS). 

We carried out point source extraction and aperture photometry of all point sources 
on the final IRAC mosaics with Phot Vis version 1.10bet a3. Phot Vis is an IDL GUI-based 
photometry visualization tool (see Gutermuth et al. 2004h that utilize s DAOPHOT modules 



ported to IDL in the IDL Astronomy User's Library (lLandsmanlll993l ). We used Phot Vis to 
visually inspect the detected sources in IRAC bands 1, 3 and 4, adding sources not detected 
automatically, but clearly visible in the images with the GUI tool and rejected any structured 
nebulosity or cosmic rays mistaken for stellar sources by the automatic detection algorithm. 
To save time, we did not visually check the band 2 photometry in this manner. Instead, we 
took the cleaned band 1 source list as the start point for finding objects in the image and 
extracted photometry at each position. Radii of the apertures and inner and outer limits of 
the sky annuli were 2.4", 2.4" and 7.2" respectively. The resultant photometry was calibrated 
using large-aperture in-flight measurements of standard stars, with an appropriate aperture 
correction in each channel to correct for the smaller apertures used in this study. In this 
paper we restrict our source catalog to only those objects with magnitude error < 0.2 in all 
four IRAC filters, a total of 18518 objects. 

We estimate the completeness in each IRAC filter by breaking up each image into a 100 
by 100 pixel grid, adding artificial stars to each grid cell and counting the number retrieved 
as a function of magnitude. Averaged over the whole W5 field, our source catalog is 90% 
complete to a magnitude of 15.5 at 3.6 /xm, 15.5 at 4.5 /zm, 14.0 at 5.8 /im and 12.7 at 8 
/im. The completeness is less in regions of bright diffuse emission: ~14 at 3.6 and 4.5 /im, 
~11 at 5.8 /im and ~9.5 at 8.0 /im. 

We conducted point source extraction and aperture photometry of point sources in the 
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Fig. 2. — W5 Spitzer IRAC bands 1, 2 and 4 color composite. Bright emission from PAH 
grains at 8 fim traces the boundary of the H ll region giving the bright pinkish color. Dense 
clusters of stars surround the O stars in the interior, with smaller clusters appearing amidst 
the diffuse PAH emission. 



Fig. 3. — Spitzer MIPS 24 /xm mosaic. The O stars labelled in Fig. 1 are marked here with 
white diamonds. Bright dust continuum emission traces the boundary of the Hll region, 
while extended emission from heated dust grains is also present in the interior (see § 3.2). 
Bright point sources are protostars. 
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Fig. 4. — Composite of emission in IRAC bands 2 and 4, plus MIPS 24 /xm. Emission at 24 
/im — presumably from heated dust — fills the interior of the Hll region cavities, see § 3.2. 
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24 nm MIPS m osaic using the point-spread-function fitting capability in IRAF DAOPHOT 
(IStetsonlll987l ). We visually inspected the image to pick out point sources not automatically 
detected due to bright diffuse emission evident throughout the image. We match the four- 
band IRAC source list to the MIPS catalog using a 2" search radius, selecting the object 
closest to the MIPS point spread function centroid in cases where more than one IRAC 
object is a match. Of the 1874 MIPS objects that match with entries in our IRAC 4-band 
list, 6 had a second object within the 2" search radius which possibly contributes to its 24 
fj,m flux. We consider this a small effect on our analysis. 



3. Analysis 
3.1. Source Classification 



In order to characterize the progress of star formation throughout W5 it is important 
to establ ish the evolutionar y status of stars within the region. The p resence of massiv e 
O stars JHillwig et all hood ) and significant molecular material in W5 (ILada et al.l Il978l ) , 
suggest that the stars associated with the region will be of relatively young age (< 10''yrs). 
Recently formed stars exhibit infrared excess emission in their spectra above that produced 
by the stellar photosphere. This excess emission arises from heated dust, either within the 
circumstellar material close to the young star, perhaps falling onto it as a part of protostellar 
collapse, or left behind after the end of star formation. This material gradually disappears 
with time, and hence the excess e naission evolves as a function of st ellar age. We make use 
of this by measuring a^R (see, e.g. iLadal 119871 : IStahler &: Pallall2005l ): 



am 



dlogX 



the value of which decreases with the progression of the star's evolutionary state, whether 
through its own aging, or the effects of its environment. 

In the case of most surveys of star forming regions, includi ng W5, infrared colors 



serve as proxies 



ample) 



o r directly measuring 



Hartmann et al 



OilR. 



Schemes developed in IWhitney et al.l (120031 ) 



Allen et al.l (l2004[). iMegeath et al.l (120041 ) and iMuzeroUe et al.l (120041 ) and tested by (for ex 



(120051) have demonstrated the power of this technique for classifying 
young stellar objects using Spitzer IRAC an d MIP S photometry. The classification scheme 
we use is described fully in lGutermuth et al.l (120081 ). We use an iterative method combining 
IRAC an d MIPS photornetry with an extinction map generated from 2MASS near-IR pho- 
tometry (jSkrutskie et al.ll2006l ). The scheme classifies the stars and filters out extra-galactic 
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contamination. 



We generate an extinction map from 2M ASS H — Kfj colors of sta rs detected across the 



2005h which itself is based 



Lombardi k Alved (120011 ) . 



entire W5 field using the method describ ed in iGuternauth et al.l ( 
on the NICE and NICER algorithms of llada et all Jl994h and 
The map has an angular resolution of ~35 arcsec, and is sensitive up to Ay ~15. However, 
the map is limited by the sensitivity of the 2MASS survey. As a result, Ay values in the map 
above ~3 may be underestimates of the true extinguishing column through the entire cloud, 
since objects behind are too faint to be detected. This is only the case for <1% of pixels in 
the map, so we ignore this effect for our analysis. With Ay values in hand for all sources, 



we der edden their IRAC magnitudes using the IR extinction law presented in lFlaherty et al. 



(120071 ). assuming each source is seen behind the full extinguishing column in each case. 



We first remove star forming ('PAH') galaxies and weak line AGN via a series of cuts in 



the fo ur band IRAC color-color diagrams after the procedure developed by lGutermuth et al. 
(120081 ). Because W5 is ~4x more distant than the nearby regions studied by Gutermuth 



et al., this filter removes many apparent young stellar objects (YSOs). In the Appendix we 
discuss how we characterize and account for this effect in our subsequent analysis. Next we 
filter out unresolved shocked blobs of P AH emission by cu tting objects with a large 4.5 /xm 
excess, i.e. very red [3.6] — [4.5] color (jSmith et al.ll2006l ). We categorize the remaining — 
presumably stellar — sample primarily relying largely on the [4.5] — [5.8] color to discriminate 
among SED classes. In Figure 5 we show IRAC color-color diagrams for sources in W5, 
marking the location of protostars (Class I/O, red dots), stars with disks (Class II, green 
dots) and stars exhibiting only photospheric colors (Class HI, black dots) as identified using 
our scheme. 

We also use 2MASS H and Ks photometry, combined with IRAC 3.6 and 4.5 /xm data 
to classify objects lacking either an IRAC [5.8] or [8.0] detection. To make sure the 2MASS 
sources have reliable photometry we require a magnitude error ^0.1 in both H and Ks 
bands. We first deredden the photometry in these four bands, and then identify IR excess 
sources as those having red [3.6] — [4.5] and Ks — [3.6] colors. The results are shown in Figure 
6, left panel. 

MIPS 24 /xm photometry provides us with additional classification information. Fig- 
ure 6 (right panel) shows the color-color diagram combining MIPS and IRAC photometry. 
Stellar sources not classified as either Class I or Class II by their IRAC or near-IR colors 
may still have red colors ([5.8] — [24] >1.5) and thus be candidate "cold disk" or "tran- 
sition disk" objects — in other wo rds: Class II objects with significant c l earing of their in- 



ner disk region (see for example , iMuzeroUe et al 



2004 



Naiita. Strom fc MuzeroUd 120071 : iFlahertv fc Muzerollell2008l : Brown et al. 



Lada et al. 



2006 



Cieza et al. 2007 



20081 ). These are 
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marked with yellow dots in Figures 5 and 6. MIPS photometry is also used as a check on the 
AG N/galaxy /shocked blob filtering, and Class I classification. AGN candidates with bright 
MIPS detection can be re-classified as protostars, and Class I objects with insufficiently red 
[5.8] - [24] or [4.5] - [24] colors are 'demoted' to Class II. 

Lastly, we also visually inspected the SEDs of all objects in Class I and Class II (a total 
of 2064 objects), where available making use of 2MASS JHKs and/or 24 /xm photometry. 
In most cases only IRAC 4-band data were present. We give the full list of objects: Class I, 
Class II, Class III and transition disk candidates in Table 2. We list source coordinates and 
photometry in near-IR JHKs (from 2MASS) and Spitzer ha.nds, and give a column denoting 
infrared SED source class and a 'ffag' column for objects with unusual SEDs. Objects are 
sorted and indexed by ascending right ascension order. Note that the full machine-readable 
table is available in the online edition of this paper. We note Class I sources with bright 
emission at 24 fim and very red IRAC— 24 /im color ([X] — [24] > 4.5, where X is the 
magnitude in any of the four IRAC bands) as a subclass: 'deeply embedded protostars' in 
the table. Table 3 gives a summary of the classification results. 

In Figure 7 we show the spatial distribution of young stars (Class I and deeply embedded 
protostars, and Class II and transition disk candidates) in W5 overlaid on the IRAC 4.5 /im 
image. Class I and Class II sources are marked in red and green respectively, transition disk 
candidates and embedded protostars with yellow and blue dots respectively. O stars are 
labelled and marked with white asterisks. 



3.2. Spatial Distribution of Young Stars in W5 

Figure 7 clearly shows that neither Class I nor Class II type objects in W5 are distributed 
uniformly, but rather (for the most part) in clustered, or filamentary formations. Within 
the larger cavity of W5W, clusters of Class II YSOs can be seen centered on the O stars HD 



17505, HD 17520 and BD +60 586 of the open cluster OCl 364 flAlter et al.lll970h . Similarly 
in the smaller W5E bubble, a cluster of young stars is visible centered on the O star HD 
18326. Smaller clusters of Class I and Class II objects are also visible in the large region of 
bright diffuse emission to the North- West of HD17505. Based on the higher proportion of 
Class I stars relative to Class H's (quantiffed in §3.4), these are probably younger clusters, 
just emerging from their parental molecular cloud. 

Extending out from the central clusters of Class II objects are more filamentary group- 
ings of stars. In W5 West, a chain of young stars appears to extend from the c omet ary bright 



rimmed clouds at the Southern rim of the bubble (object 11 in the survey of ISugitani et al. 



- 12 - 





, 1 


' ' ' 1 1 


a 


* 










■ 1 


L ■ * ' 


•al* • ■ 


[3.6] 



















-1 


■ , . . . 1 , '. . , . 





1 

[4.5] - [5.f 




-1.0 -0.5 0.5 1.0 1.5 2.0 
[5.8] - 18.0] 



Fig. 5.— Left panel: [3.6] - [4.5] vs. [4.5] - [5.8]; right panel: [4.5] - [5.8] vs. [5.8] - [8.0] 
IRAC color-color diagrams used for identifying candidate protostars. Black dots — Class III, 
green — Class II, red — Class I, yellow — transition disk candidates. AGN and PAH galaxy 
candidates are not included here. 
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Fig. 6. — Left panel: dereddened 2MASS Ks vs. Ks — [4.5] color-ma gnitude diagram. 



Dash ed red box shows location of AGN detected by IRAC Shallow survey (lEisenhardt et al. 



20041 ). Right panel: IRAC and MIPS [3.6] - [5.8] vs. [8.0] - [24] color-color diagram. 
Color scheme: black dots — Class III, green — Class II, red — Class I, yellow — transition disk 
candidates. Suspected AGN and PAH galaxies are not included here. 



Table 2. W5 Source List 





RA 


Dec 


7 

J 






[■^ fil 


[4.0J 


si 


[».UJ 


[Z4J 




ID 


(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Type'' Flag'' 


1 


41.081526 


60.510607 


15.34(06) 


13.69(03) 


13.04(03) 


12.60(01) 


12.50(01) 


12.36(04) 


12.47(07) 




III 


2 


41.098560 


60.682772 


10.90(02) 


10.68(02) 


10.58(02) 


10.50(01) 


10.49(01) 


10.42(01) 


9.96(02) 




III 


3 


41.100737 


60.667424 


13.30(02) 


12.26(02) 


11.88(02) 


11.63(01) 


11.72(01) 


11.43(04) 


10.90(11) 




III 


4 


41.109056 


60.660875 


13.01(02) 


12.52(04) 


12.29(03) 


12.18(01) 


12.12(01) 


11.81(05) 


10.67(05) 




II 


5 


41.109473 


60.515686 


12.52(02) 


12.09(02) 


11.98(02) 


11.94(01) 


11.94(01) 


11.88(02) 


12.02(05) 




III 


6 


41.112052 


60.519115 


13.27(02) 


12.53(02) 


12.12(02) 


11.95(01) 


11.86(01) 


11.75(02) 


11.74(05) 




III 


7 


41.116165 


60.509403 


14.20(03) 


13.73(03) 


13.53(04) 


13.48(01) 


13.52(01) 


13.54(08) 


13.43(17) 




III 


8 


41.128907 


60.512183 


14.90(04) 


14.12(04) 


13.67(05) 


13.41(01) 


13.32(01) 


13.36(08) 


13.78(19) 




III 


9 


41.128946 


60.536046 


13.96(03) 


12.39(02) 


11.74(02) 


11.36(01) 


11.35(01) 


11.20(01) 


11.12(03) 




III 


10 


41.130516 


60.673029 


12.78(03) 


12.30(03) 


12.03(03) 


11.89(01) 


11.82(01) 


11.74(02) 


11.70(05) 




III 



Note. — Table 2 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. Values in parentheses by photometry signify error in last 2 digits of magnitude value. Right ascension and Declination 

coordinates arc J2000.0. 

"Source class as defined in text. 

''Sources with unusual SED shape flagged with 'v' or uncertain source class '?' 
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Fig. 7. — Source distribution overlaid on the MIPS 24/im greyscale image. Color scheme: 
red — Class I, green — Class II, yellow dots — transition disk candidates. O stars marked with 
asterisks. 



- 15 - 



199ll ) up to the clusters around HD 17505 and HD 17520, and through the cluster around BD 
+60 586 to the bright rimmed clouds on the Eastern rim of the bubble. In W5 East, a chain of 
youii g stars extends East- West from each of the two peaks of AFGL4029 (jPrice &: Murdoch 
19831 ) through the cluster around HD 18326 to the bright PAH ridge that marks the Western 
edge of the bubble. 

Numerous protostars (Class I objects) are seen in projection on the cloud rim at the 
PAH emission/Hll region cavity boundaries. For example, the bright rimmed clouds that 
make up the eastern border of W5 East, the bright rimmed clouds at the North-Eastern rim 
of W5 West, and the cometary clouds at the southern rim of W5 West all contain several 
Class I objects. Several isolated cloud remnants and small 'elephant trunk' formations within 
the Hll region cavity also contain Class I or Class II objects — these are listed in Table 4. 
We give a 'type' to each object to distinguish the different morphologies. These objects are 
interesting candidates for triggered star formation on small scales, see the discussion in § 
4.3. Figure 8 shows several examples of these objects. 

The influence of star formation on its environment is clear to see in both the IRAC 
(Figur e 2) and MIPS irn ages (Figures 3 and 4). Bright diffuse emission in the four IRAC 
filters (jPraine fc Lil 120071 ) is only present around the edge of the two main Hll region bub- 
bles in roughly ring-like structures, and brightly lit-up over the molecular clouds W5NE 
and W5NW. Presumably the PAHs responsible for much of the emission in IRAC bands 
1, 3 and 4 are destroyed in the harsh ionizing environment of the Hll region. At 24 fim 
the morphology is similar in the bubbles' rims, but very different diffuse emission is seen 
in the interior surrounding t he 5 O stars (as labelled in Figure 3). See n in n umerous mas- 
sive star forrn i ng re gions by IChurchwell et al.l (120061 ) , ISmith fc Brooks! (120071 ) and recently 
Harvey et al.l (120081). this ad ditional component of diffuse emission is consistent with the 
picture of ICrete et al.l (119991 ) whereby emission within Hll regions at mid-IR wavelengths 
is produced by heated larger dust grains (radius > 12 A) that survive being destroyed by 



Table 3. Source classification summary 



Class 


N Objects 


Class I/Protostars 


176 


Class II 


1806 


Class Ill/Photospheres 


15709 


Trans. Disk 


79 


Embedded Protostars 


3 


AGN 


729 


PAH Galaxies 


198 


Total 


18518 



Fig. 8. — Examples of star forming globules in W5 shown at 8 fim. PAH emission is seen at 
or near the globule heads. Objects are labelled as in Table 4. 



ionizing radiation. As iDraine &: Lil (120071 ) have shown, emission in IRAC bands 1 and 4 is 
much more sensitive to the PAH mass fraction than is emission in the MIPS 24/im band. 
Further, 24 /xm band emission increases rapidly when the interstellar radiation field (ISRF) 
rises to 10^-10^ times the local value. Thus a region with low PAH mass fraction but high 
ISRF can have enhanced 24 /im emission, but will have low IRAC 3.6 and 8.0 fim PAH 
emission. The morphology of the 24 /xm emission suggests that these larger grains in the 
immediate vicinity of the O stars are blown away by the strong stellar winds. 



3.3. Spatial Variation in SED types 

Class I and II objects are present in varying amounts across W5. Although with our 
present dataset we cannot tell if any one Class I object is younger than a given Class II 
object, these classes do represent different evolutionary stages in the formation of stars. In 
W5, we find that in regions where our ^^CO emission map exceeds T^ = 7.5 K km s~^, there 
are on average 0.23 Class I objects for every Class II. In regions with T^ < 7.5 K km s^^, 
the ratio is 0.04. This suggests that the stellar clusters associated with significant molecular 
material in W5 are systematically younger than those in cleared out regions, for example 
the Hll region cavities. 

The varying levels of completeness in each IRAC map as a function of association with 
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the bright, diffuse emission affect our result. As described in § 2, the 90% hmiting magnitude 
decreases by ~1.5 mag in IRAC channels 1 and 2, by ~2.5 mag in channel 3 and by ~3.2 
mag in channel 4 in going from the e vacua ted H ll region to the brightest regions of diffuse 
emission. Recently, IChavarria et al.l (120081 ) studied the effects of reduced sensitivity in the 
5.8 and 8 fim Spitzer bands, combined with the obscuring effect of bright, diffuse PAH 
emission on the detection of Class I and Class II sources in IRAC surveys. They found that 
Class II sources are on average fainter than Class I's, and are likely preferentially missed in 
bright PAH emission regions as a result. 



3.4. Clustered versus Distributed Star Formation in W5 

We want to describe in an objective way how the stars in W5 are grouped and address the 
following question: is there a population of stars that do not belong to identifiable groups 
that may have formed in isolation, in other words a 'distributed' population? Lacking a 
dataset that describes the masses and dynamics of the stars and gas which would allow us 
to assign stars to clusters based on gravitational association, we base cluster membership 
only on spatial arrangement on the sky and so we do not require that these clusters be 
gravitationally bound. We restrict ourselves to infrared excess sources only to simplify the 
issue of whether stars are associated with W5. We include the deeply embedded sources 
and transition disk candidates along with the Class I and Class II objects and exclude stars 
exhibiting only photospheric emission (Class HI sources), AGN and PAH galaxies for a total 
of 2064 objects. The source distribution is shown in Figure 7, and reproduced in Figure 9, 
without the background image. 

Determining a criterion for cluster membership from spatial distributions alone is some- 
what arbitrary. In this paper we use the so-called 'minimal span ning tree' (MST) method 



to id e ntify and characteri ze clustering in W5 (see recent work by ICartwright &: Whitworth 



2004J : iBastian et al.l 120071 ). iGower fc Rossi (Il969l ) describe how such a tree is constructed 



from a list of source positions. To construct the tree, we first generate the network of lines 
that joins together the positions of all objects in our input list, such that the total length of 
all lines joining points is minimized and there are no closed loops. Each object is assigned 
one 'branch length,' that is: the projected distance to its nearest neighbor. We plot the 
distribution of branch lengths in Figure 10. 

To distinguish 'clustered' objects from 'distributed' using the MST requires 2 things: 
a break/cutoff length, ds and a minimum group size A^. A group is then a collection of 
objects linked by branches shorter than dg with at least stars (we choose A^ = 10). We 
can estimate a value for dg from the distribution of branch lengths in Figure 10. Following 
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Fig. 9. — Distribution of stellar sources as in Figure 7. Green points: Class II, red: Class I, 
yellow: transition disk candidates. Black asterisks mark the location of the known O stars. 
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Object Number 

Fig. 10. — Minimal spanning tree length distribution, showing each object's branch length 
(y-axis) versus its number in the sorted list. Straight line fits (dashed lines) to the short and 
long length parts of the distribution intersect at 0.86 pc. 



Gutermuth et al. (2008, in prep) we fit straight lines through the long and short-length 
portions of the distribution. Where these lines cross defines the break length for clustering. 
For the IR excess sources in W5 this method yields dg = 0.86 pc. Using this break length 
we find 16 groups containing at least 10 stars and a clustered fraction of 70%. 

A single value of dg chosen in this way may not be representative of the varying spatial 
density of stars across the region, and may underestimate the amount of substructure within 
groups. Altering dg changes the number of groups and the relative fraction of clustered 
and distributed sources. As the cutoff length is increased, an increasing fraction of objects 
qualify as 'clustered.' The number of groups (Ngrp) we find also increases up to a maximum 
before falling off as groups start to merge until all groups belong to a single 'cluster.' These 
trends are shown in Figure 11. An alternate value for the break length (i,, can be derive d 



from the peak in the plot of Ng^p versus increasing dg. As argued by iBattinelli I (Il99ll ). 
this method can be thought of as returning a maximum of information from our source 
distribution. In this case we find a cutoff of dg = 0.54 pc, Ngrp = 27, and a clustered fraction 
of 44.2%. In Figure 12 we compare the groupings produced by the two methods side by side 
to demonstrate this. 



Both values of dg used here ar e large when cornpared to the typical separations of stars 
in nearby regions such as Taurus (IHartmann et al.l 120051 ). where a cutoff length of 0.54 pc 
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would incorporate many objects not associated with any cluster. However, Taurus is a low- 
mass star forming region and so the characteristic scales of clustering are likely to be smaller 
there. In recent surveys of low mass s tar forming region s (lEnoch et al.l 120071 ) and the high 
mass star forming regions in Cvgnus (IMotte et al.l 120071 ). typical star- forming dense cores 
are ~0.1 pc in size. However, jMotte et al.ll2007h also found that massive dense clumps can 
be much larger: 0.5-0.8 pc, which may be a more relevant scale for cluster and massive star 
formation than for low density, more isolated star formation seen in Taurus. 

Table 5 summarizes the clustering results from the two methods. The value for the frac 



tion o f stars in clusters from both these methods is lower th an that reported by lAUen et al. 



( 120071 ) in Spitzer sniyeys of Orion and Ophiuchus (74-78%). ICarpenterl (120001 ) . in his analy- 
sis of the Orion A and B, Mon R2 and Perseus molecular clouds with 2MASS found a range 
of clustering fractions, different for every cloud, between 56-100%. This fraction could be as 
much as a factor two lower depending on how the age of the distributed population is mod- 
eled. In W5 — at a distance of 2 kpc — we are sampling much less of the IMF than in these 
nearby regions, which lie within 1 kpc of the Sun. We also do not know the disk-fraction of 
YSOs in the W5 clusters. ICieza et al.l (120071 ) have shown that a significant fraction of stars 
lose detectable circumstellar disks (in the mid-IR) on a timescale ~1 Myr. These objects 
would not be included in our sample of IR-excess sources. 

The mass completeness limit for our survey is dictated by our photometric completeness 
and the wide span of YSO mid-IR colors. The primary constraint in completeness for the 
IRAC 4-band sample is the 8 /xm detection limit, since this has low sensitivity and very 
bright background contamination from PAH feature emission. A typical photosphere with 
[3.6] — [8.0] ~0.1 at age 2 Myr, and [8.0] = 12.7, equates to a limiting mass ~2 Mq (using 
the evolutionary models of D'Antona & Mazzitelli 1994, and adopting Ay = 2, a typical 
value for the Hll region cavity: Hillwig et al. 2006). On the bright diffuse emission the 
8 fim detection limit for a photosphere equates to ~8 Mq. The additional 2MASS+IRAC 
bands 1 & 2 sample is limited by the Ks band comp leteness. We e stima te that the 90% 
completeness in this band is 14.2. Using the models of iBaraffe et al.l (Il998l ). this equates to 
a star of 1 Mq at 2-2.5 Myr, at an Ay of 2 and a distance of 2 kpc. Both the 8 /xm and Ks 
band mass completeness limits are conservative estimates, since YSOs with infrared excess 
emission will be detected more readily. For example, the [3.6] — [8.0] color peaks at ~1.4 for 
IR excess YSOs in W5. In this case, [8.0] = 12.7 converts to a mass of ~0.8 Mq at 2 Myr 
in the cavity. In summary, we estimate that for a young population of 2 Myr age, we are 
complete to ~8 Mq for photospheres and ~4 Mq for typical disk excess objects in the IRAC 
4-band sample seen against the bright background. In the cavity we are complete to ~2 Mq 
for photospheres and ~0.8 Mq for disk excess objects. In the additional 2MASS sample we 
are complete to a mass ~1 Mq. 
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Table 4. Star forming globules in W5 



Object RA (deg) Dec. (deg) Type 



1 


2/t 


49"" 55? 78 


60° 


26' 


24^' 46 


comet 


2 


2h 


SC" 28? 16 


60° 


09' 


30^44 


comet 


3 




50™ 33? 85 


60° 


25' 


47'.' 12 


comet 


4 


2h 


51™ 32?90 


60° 


03' 


53'.' 34 


trunk 


5 


2h 


51™ 53?91 


60° 


06' 


56'.' 97 


comet 


6 


2h 


52m 17s 19 


60° 


03' 


17'.'22 


trunk 


7 


2h 


52m 2if75 


60° 


54' 


11'.' 86 


trunk 


8 


2h 


53™ 39? 14 


60° 


46' 


50^'51 


trunk 


9 


2h 


53™ 39?22 


60° 


30' 


17" 95 


comet 


10 


2h 


54™ 52? 34 


60° 


35' 


43'.' 16 


trunk 


11 


2^ 


54™ 58? 35 


60° 


41' 


43'.' 76 


trunk 


12 


2h 


56™ 8? 69 


60° 


10' 


25'.' 26 


trunk 


13 


2h 


58™ 37? 57 


60° 


41' 


53'.' 69 


comet 


14 


2h 


59™ 46? 22 


60° 


21' 


10'.' 32 


comet 


15 


3'' 


0™ 23? 53 


60° 


17' 


55^' 05 


trunk 


16 


3*^ 


1™ 1?98 


60° 


21' 


57('54 


trunk 


17 


3'^ 


1™ 17?47 


60° 


24' 


13" 20 


trunk 


18 


■-)'' 


1"' 17^91 


()()■•' 




2l'.'():-! 


trunk 



Note. — Coordinates are J2000.0. 



Table 5. Clusters/Groups in W5 



Method 

Straight-line fit Ngrp maximum 



Total Ni, in clusters 


1444 


912 


Group size (iV*) 


10-459 


10-201 


% in clusters 


70.0 


44.2 


Ni/Nii 


0-50% 


0-50% 



-23- 




Fig. 11. — Minimal spanning tree results. Solid line and left y-scale: clustered fraction of 
stars (belonging to all groups >10 stars) as a function of branch length cutoff. Dashed hne 
and right y-scale: number of groups {Ng^p) containing 10 or more stars identified by the MST 
algorithm. NOTE: coarse steps in x-axis miss finer detail in Ng^p] true peak is 27 groups at 
ds = 0.54 pc. 
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Fig. 12. — Upper panel: the MST groups as found using the straight-hne fit technique: ds 
= 0.86 pc, Ngrp = 16. Stars are colored according to their group for visual identification 
only. Lower panel: MST groups picked out at the maximum in Ng^p vs. dg plot, using dg = 
0.54 pc, Ngrp = 27. In both plots, black dots represent stars not associated with any cluster. 
Black asterisks mark the location of the known O stars. 
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Discussion 



4.1. Origin of the Distributed Population 

Our clustering analysis shows that there is a significant distributed population of young 
stars (30-57% of IR excess objects). These may have formed in groups or clusters and 
dispersed through random motions or been ejected at high velocity. Conversely, they may 
have formed in their current locations in largely unrelated, isolated events. 

On a ~10 Myr timescale, the distributed population may be a natural consequence of 
dynamical interactions between young stars in groups throughout the original giant molecular 
cloud. Small groups {N < 36) have short relaxation times relative to their crossing times 
and, once the gas is re moved, will not remain visible as clusters for much longer than this 
( lAdams fc Myerd l200ll ). However, if we assume that the distributed stars are only as old 
as the oldest O star in W5 (a reasonable assumption, given that we have identified them 
through their IR excess emission), they will have had < 5 Myr to travel. 

Extreme dynamical interactions, namely close encounters between stars in multiple 
systems, can produce high velocity runaway stars. Numerical studies show th at a given 
triple s ystem will eject a membe r within about 100 crossing times (~30000 years: iReipurth 



2000l ). ISterzik fc DurisenI (119951 ) calculate similar timescales for close systems consisting of 
5 stars, and predict velo cities of 3-4 km s~^ for the ejected stars. Few have been found; 
Goodman fc Arcd (120041 ) find 7 high velocity (>10 km s~^) candidates in the literature. 



Fiiresz et al. (2006, 2008) measured the radial velocities of young stars in NGC 2264 
and the Orion Nebula Cluster, and found dispersions of 3.5 and 3.1 km s~^, respectively. In 
W5, the average projected distance of distributed stars from the nearest cluster is 8.5 pc. 
Assuming a velocity of 3 km s^^ perpendicular to the line of sight, this distance could be 
traversed in ~3 Myr — well within the current upper age upper limit for W5 (5 Myr, see: 



Karr fc Martinll2003l ). Thus it is plausible that at least some of the distributed population 
originated in nearby clusters. 



4.2. Star Formation Efficiency 

We calculate the star formation efficiency in W5 using the equation: 



stars 



Mstars + Mgas 



(2) 
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Here Mgtars represents the mass in young stars only, as identified by their infrared excess, so 
e is a current efficiency, averaged over the last few million years. For the entire W5 region we 
count 2064 young stars and assume an average mass of 1 Mq. To find the gas mass we use 
the Ay map of W5 described in § 3.1 and estimate the molecular gas colu mn density from th e 



standard relation: N{H2) = 1.9 xlO'^^Ay protons cm ^ (see, for example: iBohlin et al.lll978l ). 



Assuming all gas is associated with W5, we convert N{H2) to a project ed two-dimensional 



mass density at each pixel in the Ay map: k = 15 xAy Mq pc ^ (as in iLombardi fc Alves 



200ll ). Each Ay map pixel has a size 0.34x0.34 pc. Over our whole survey field, we find a 



total gas mass of 6.5x10^ Mq and derive an efficiency of 3%. 

We calculate e for the clusters identified by our MST treatment in § 3.4. For each of the 
16 clusters, we count the number of stars contained, assume a stellar mass of 1 Mq and add 
up the gas in the corresponding pixels. Inside the H ll region cavity, where the molecular gas 
has been destroyed, the efficiencies are high (26-39%). In the rim of W5, where molecular 
gas remains and embedded clusters are found, e is 10-17%. These latter values are likely 
a lower limit, as our census of young stars is incomplete in the bright extended emission 
coincident with the embedded clusters. 

Recent results from the c2d survey (Evans et al. 2008 in prep.) show that nearby dark 
clouds where low mass star formation is occurring have star formation efficiencies of 2-4% 
cloud-wide and 15-20% in clusters, similar to the values we have obtained for the high-mass 
star forming reg i on W 5. Our results for e in the clusters are also comparable to those of 



Chavarria et al.l (120081 ) for the S255 region 



4.3. Triggered Star Formation in W5 



We consider two mechanisms of triggered star formation relevant to W5. The first 
mechanism, called radiatively driven implosion (RDI), is the compression of pre-existing 
density enhancements (small cloud clumps or globules) inside and on t he boundary of th e 
ionized bubble by the high pressure of the Hll region (for a review, see iKlein et al.lll985l ). 
Clumps of material visible as bright PAH emission in IRAC bands 1, 3 and 4 with spatially 
coincident Class I or Class II objects are summarized in Table 4. These include isolated 
clumps, 'el e phant trunk' formations and bright rims (as shown in Figure 8). As discussed in 
Elmegreenl (119981 ) . if globule squeezing does occur, its timescale is expected to be short, since 
an isolated over density can be co r apress ed immediately on being engulfed by the H ll region 
or stellar wind. [Thompson et al.l (120041 ) investigated the cometary features on the southern 
rim of W5 West for signs of star formation triggered via the RDI mechanism (globules 4-6 
in our Table 4). They found that the Ha, CO-molecular and dust morphologies of the three 
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pillars are 'reasonably consistent' with the model of iLefloch fc Lazareffl (119941 ) for radiatively 
driven implosion, at the early collapse phase. They used near-IR photometry to detect several 
candidate young stellar objects — these we confirm through our Spitzer photometry as Class 
I and Class II objects. They estimate that the timescales for duration of UV illumination, 
shock crossing ti r nes ac ross the features, and protostar/YSO ages are all ~10^ yrs. Since 



Thompson et al.l (120041 ) were only able to establish that approximate pressure equilibrium 
holds in the pillar heads from their data, we need a more detailed study of the gas dynamics 
in these pillars to determine if they are currently collapsing due to external pressure from 
the Hll region. 

The second mechanism, 'collect and collapse,' is the large scale expansion of an ionized 
bubble, powered by an over-pressurized Hll region and stellar winds. This expansion can 
drive shock fronts into the surrounding medium, sweeping up a d ense ridge or shell which 
collapses into stars when a critical density of material accumulates. IWhitworth et al.l (119941 ) 
presented an analytical treatment of this process. A shock front forms and gathers material 
until it is able to fragment and collapse to form stars, at a time t/rag and radius Rfrag given 
by: 



tfrag = 1.56 Myr aj^^ L 



1/11 -5/11 



49 



n 



Rfrag — 5.8 pC Oq 3 L 



4/11 ^1/11 -6/11 



49 



Tin 



(3) 
(4) 



where ao.2 is the sound speed inside the shocked layer in units of 0.2 km s~^, L49 is 
the central source ionizing flux in units of 10^^ photons s~^ and ns is the initial gas atomic 
number density in units of 10^ cm~'^. 

In any massive star forming region, presumably some combination of H II region expan- 
sion and stellar wind forces operate. Stellar winds certainly play a significant r ole in shaping 
the rn orph ology of an H I I regio n only until an age of ~ lO^years. As shown by lMcKee et al. 



(119841 ) and lWeaver et al.l (Il977l ). weak stellar winds (wind luminosity <C 1.26x10^^ x L49 x n, 
where n is the ambient density) are likely to be confined by the Hll region pressure to a 
smaller bubble around the ionizing source, since gas that is evaporated from clumps of ma- 
terial in the H ll region mixes with the hot stellar wind bubble and can radiate energy away 
efficiently. Strong winds can overcome this confinement — they create a bubble with size de- 
termined by the photoevaporation and displacement of surrounding inhomogeneities in the 
ambient gas. This then expands along with the Hll region. In estimating the lifetimes of 
the W5 bubbles, we thus consider only the effects of the H ll region expansion. 

As described in §1, the gross morphology of W5 in the mid-IR is defined by two large, 
roughly circular rings of bright PAH emission that mark the smaller W5 East (W5E, radius 
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~10 pc), and larger, more irregular W5 West (W5W, overall radius ~22 pc). At the approx- 
imate center of W5E lies HD 18326, an 07V star surrounded by a dense cluster of young 
stars. It seems reasonable to assume that whether this is a sphere or a ring of gas, that the 
cluster lies near the true center, and is not si gnificantly offs e t eith er to the foreground or 
background. Adopting stellar parameters from lMartins et al.l (120051 ). an 07V star will have 
ionizing flux L49 = 0.43. Since we cannot calculate appropriate values for ao.2 and ns from 
our current datasets, we adopt 1.0 for b oth. The H ll region fra gmen tation time we calc ulate 
is 1.69 Myr at a radius of 5.37 pc. Both lwhitworth et al.l Jl994l ) and bale et aP J2007ah note 
that the density in a real region is likely to be lower — this would make the fragmentation 
time longer and the radius larger by a factor ~2. 

W5W presents a less clear-cut picture. It contains an isolated 08V star: HD 237019, 
as well as three dense clusters centered on BD +60 586 (07.5V), and the multiple systems 
HD 1750 5 (two 07.5V(('f)) s tars, an 06.5III((f)) and an 08.5V) and HD 17520 (09V + 
Be). See iHillwig et al.l (120061 ) for a more detailed study of the O stars in W5. The relative 
configuration of these four systems, and their relationship to the surrounding diffuse gas 
(whether in front or behind) is not known. We consider here two simple scenarios and 
assume that all four objects and the ring of PAH emission are in the same plane. Scenario 1: 
HD 237019 represents an initial episode of star formation, and its isolation is due to its cluster 
having had time to disperse. In the Whitworth model, an 08V star has a fragmentation 
time of 1.8 Myr and collapse of swept up gas occurs at a radius of 5.0 pc. In this picture, this 
event triggered the formation of the three dense clusters in W5W which all lie at roughly 
6.8 pc from HD 237019. These clusters, and the continuing expansion of the Hll region 
then triggered the current, ongoing star formation seen in association with molecular clouds 
W5NW and W5NE. Scenario 2: the three dense clusters form together — with their present 
arc-like arrangement due to some initial filamentary distribution of molecular material. They 
then trigger star formation in the remaining molecular material to the north, north-east and 
south. The 07.5V star BD +60 586 by itself has a fragmentation time of 1.75 Myr and 
collapse occurs at a radius of 5.16 pc. The corn bined systems HD 17505 and HD 17520 
(with at least 5 O stars, see: iHillwig et al.l 120061 ) have a fragmentation time of 1.46 Myr 
and collapse occurs at a radius of 6.19 pc. In this picture, HD 237019 would have been 
ejected from one of the multiple O star systems. Located at ^10 pc from the nearest cluster, 
assuming a velocity entirely transverse to our line of sight of 10 km s~^, it would take ~1 Myr 
to arrive at its curre nt location, consistent with the likely young age of the cluster as a whole. 
Moffat et al.l (119981 ) found a runaway frequency percentage of 14% among Galactic O stars 
and included this star as a marginal runaway ca ndidate. Alth ough the most favored creation 
mechanism for runaways is supernova ejection (jBlaauwlll96ll ). ejection from compact young 
stellar clusters or through binary interaction is also thought to be responsible for some O 
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star runaways flPoveda et all 119671 : IClarke fc Pringldll992f ). Moffat et all fll998[ ) found HD 
237019 to be only a marginal runaway candidate, given the large uncertainty in its tangential 
velocity: 15±17 pc Myr~^. 

There are several observations consistent with the collect and collapse mechanism in W5. 
The low lev el of non-thermal r adio emission suggests there have been no supernovae during 
its lifetime (IVallee et al.lll979l). The m ain sequence lifetimes of the O stars are estimated 
to be < 5 Mvr jKarr fc Marti n"2003V The timescales for triggering presented above are 
shorter than this. IWilking et al.i (il984i ) found several young embedded OB stars around W5 
in the molecular clouds, and argued that their locations near to the cloud edge, together 
with some evidence of magnetic field alignment p arallel to the i o nizati on front are suggestive 
of triggering via external compression. Finally, iNakano et al.l (120081 ) found evidence of an 
age difference of ~3 Myr between the cluster around HD 18326 in the W5 East bubble and 
the young stars in the rim. This is longer than our estimate above of ~1 Myr, but given 
the uncertainties in gas sound speed a and initial cloud number density n, is an allowable 
timescale in the collect and collapse model. 



Karr fc Martini (120031 ) argue that the distribution of young stars in W5 indicates a 



shorter timescale, more consistent with RDI than with collect and collapse. They investigated 
the locations of IRAS point sources with colors corresponding to YSOs. Their results did 
suggest a 2-generation model, with older stars at the centers of the H ll regions and ongoing 
star formation surrounding them, however they found that the number density of young 
objects peaks ~5 pc inside the Hll region (as defined by the 6.2K 1420 MHz contour), which 
corresponds to a triggering timescale of 0.5-1 Myr. Of the 42 IRAS sources from their list 
that fall within our Spitzer image we find only 16 in our source list at 24 /im. The remainder 
may be spurious (they are also not seen in our 70 fim image) and in fact may be small knots 
of IR emission that masquerade as point sources in the IRAS survey. All are associated 
with bright PAH emission, includi ng one coincident wit h globule 8 in Table 4. The peak 
in IRAS-classified YSOs found by iKarr &: Martini (120031 ) is also not upheld by an analysis 
of our Spitzer sample which shows that the Class I objects are found predominantly along 
the cloud rims or at the cloud centers. This distribution, with the YSOs further from the O 
stars argues for the longer timescale of collect and collapse. 

Studies of a young embedded object G138.295+1.5 55 and the UGH ll region G138.300+1. 558 



in AF GL4029 (W5E) with near- and mid-IR imaging (jPeharveng et al.lll997l : IZavagno et al. 



19991 ) suggest that collect and collapse is not playing a role in star formation here, since 
the youngest object (the former) is apparently closer to the ionization front than the older, 
latter object. However the true 3-dimensional configuration of all the objects involved is not 
exactly known, this may be due to projection effects. The clustering of Glass II objects im- 
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mediately outside the cloud to the west of G138. 295+1. 555 and G138. 300+1. 558 is certainly 
suggestive that the sequence of star formation here is west to east, i.e. outside-in, away from 
the ionizing star HD 18326. 

In consideration of our simple model and the results in the literature, it seems plausible 
that both the RDI and collect and collapse mechanisms are at work in W5; RDI on the smaller 

scale of comctary globules, and collect and collapse on the larger scale of the Hll region. 
Detailed investigations (Kocnig et al. 2008, in preparation) of the spatial distributions and 
relative ages of the YSOs may further constrain the scenarios presented here. 

5. Conclusions and Future Work 

We have presented initial results from our extensive Spitzer survey of W5. Shorter 
wavelength data from IRAC (3-8 //m) and longer (24 fxia) wavelength data from the MIPS 
instrument were combined to maximize spectral coverage of detected sources. 

Even before photometric analysis, dense clusters of stars are clearly visible across the 
region, centered on the O stars HD 18326, BD +60 586, HD 17505 and HD 17520, and also 
across the extensive PAH emission that defines the outline of W5. At 24 /xm, substantial 
extended emission is visible, presumably from heated dust that survives in the strongly 
ionizing environment of the H II region. 

We used photometry of more than 18000 point sources to analyze the spatial distribu- 
tions of young stars, establish their evolutionary status via their infrared colors and mag- 
nitudes, and assess their clustering properties across this large star forming region. The 
large clusters that dominate the region, centered on the massive O stars, contain numerous 
infrared excess sources. Looking at the large scale distributions of stars at different evolu- 
tionary stages, we find that within the evacuated cavity of the H ll regions that make up W5, 
the ratio of Class II (older) to Class I (younger) sources is ~7 times higher than for objects 
detected coincident with the molecular clouds in the rim. We attribute this difference to 
an age difference between the two locations, and consequently postulate that at least two, 
distinct generations of star formation are visible in the region. An isolated O star in W5 
West, HD 237019, may represent an initial episode of star formation in the region, preceding 
the formation of the large clusters in the cavity, although we cannot rule out its ejection 
from an O star multiple system, for example HD 17505. 

The clustering results show that, considering infrared excess sources alone, (2064 ob- 
jects) ~45-70% are found in clusters with >10 members. Incorporating the sources appar- 
ently misclassified as AGN and PAH galaxies (see Appendix) extends this range to ~40-70%. 
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The remainder are in the distributed mode, many of which could have formed in nearby clus- 
ters. 

We looked at the role that triggered star formation may have played in W5. We cat- 
alogued isolated globules of diffuse PAH emission, and so-called 'elephant-trunk' structures 
that contain young protostars or stars with protostellar disks. These are examples of pos- 
sible RDI-trigg ered star formation eve nts. On the larger scale, we tested the analjd;ical 



formulations of IWhitworth et al.l (119941 ) for the collect and collapse mechanism. Our sim- 
ple estimates show that triggering remains a plausible mechanism to explain the multiple 
generations of star formation in W5 and merits further investigation. 

Substantial work is ongoing to refine our understanding of star formation in W5 in the 
light of triggering models and the clustering of stars. We have undertaken a near infrared 
survey of the whole region in J, H and Ks bands to extend the stellar SEDs to shorter 
wavelengths and detect more young stars against the bright background emission in W5. 
This will permit a full analysis of the clustering of young stars in the region. We have 
obtained optical spectra of several hundred stars across W5 with the aim of determining 
their relative ages by constructing H-R diagrams. This will allow us to better understand the 
history of star formation across the entire region. A much clearer picture of the progression 
of star formation across W5 and the role of feedback in this process should result. 
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Appendix 
A. AGN/PAH Galaxy filtering 



Following iGutermuth et al.l (120081 ) PAH galaxy candidates are identified on the basis of 



the following criteria: 

1 0^ 

[4.5] - [5.8] < ^([5.8] -[8.0]-!), 

[4.5] - [5.8] < 1.05, 
[5.8] - [8.0] > 1, 

[3.6] - [5.8] < ^([4.5] - [8.0] - 1), 
[3.6] - [5.8] < 1.5, 
[4.5] - [8.0] > 1. 

After removing these sources, AGN are then picked out following these criteria: 

[4.5] - [5.8] > 0.5, 

[4.5] > 13.5 + ([4.5] - [8.0] - 2.3)/0.4, 
[4.5] > 13.5. 

We plot the distribution of 729 AGN and 198 PAH galaxy candidates in W5 in Figure 
13. It is notable that these objects show obvious non-uniformity. The non-uniformity arises 
from several sources: 1) variation in the intrinsic distribution of galaxies toward (behind) 
W5, 2) variable extinction across the region due to the molecular clouds that make up W5 
and material along the line of sight through the Galaxy in this direction, 3) variable levels 
of photometric completeness across the field due to bright extended nebular emission and 
4) clustering coincident with the positions of the dense clusters of young stars around the 
O stars in the W5 Hll region cavity, due to faint (low mass) stars misclassified as AGN 
or PAH galaxies. This last issue affects all of our analysis, thus we need to characterize — 
statistically — the properties of these misclassified objects and how adding them back to the 
original sample affects our previous results. 

We consider the AGN first. In the extreme case where all objects are misclassified, we 
return the full list of 729 AGN to the stellar list, classify the returned objects into Class I, 
II etc., and examine the clustering fraction. With the MST straight line fit method, we find 
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ds = 0.88 pc, N, 



grp 



Ngrp distribution we find d, 



21 and a clustered fraction = 68.6%. If we use the maximum in the 
0.5 pc, Ngrp = 26 and a clustered fraction = 40.4%. 



Instead of this, presumably some intermediate fraction of AGN should be returned as 
AGN. We must first calculate the expec ted number and a veraged distribution of AGN in 
W5, using a sample of AG N extracted by IStern et al.l (120051 ) from the Spitzer IRAC shallow 
survey of the Bootes field (lEisenhardt et al.ll2004l ). convolved with our extinction map (as 
described in § 2.1) and with completeness estimates derived for our Spitzer survey of W5 via 
a simple Monte Carlo simulation c ode. T he model AGN sample is drawn from a 7.7 square 
degree section of the IStern et al.l (120051 ) survey, identifyin g sources as AGN base d on the 
same color classification scheme that we use in this paper (iGutermuth et al.ll2008l ). 



We use the IRAC completeness maps as described in § 3.1. The area imaged in our 
W5 survey (3.5 square degrees) is smaller than that covered by the Bootes data. In our 
Monte Carlo simulation, we input a proportionately smaller, random selection of AGN from 
the Bootes list and distribute these randomly in the field of W5. Each object is assigned 
an Ay value based on its location in the extinction map. We appl y this extinctio n to it s 
IRAC photometry using the infrared extinction relation presented in iFlaherty et al.l (120071 ). 
The corresponding location within the completeness maps gives us the completeness value 
(a number <1) at each of the four IRAC wavelengths given its extincted magnitude. A 
random number generator is then used to determine whether or not an object is detected, 
with the completeness estimates as a measure of the probability of detection. For detection 
we require that an object be detected in all four IRAC bands. This simulation predicts on 
average 270 AGN detected in W5. 

This distribution of AGN has a certain characteristic space den sity distribution. At eac h 
iteration we measure the 'nearest-neighbor' distance distribution (ICasertano fc Hutlll985l ). 
For each object we find the projected distance to its sixth nearest neighbor in arcsec (dg) and 
construct a histogram of distances. The final output of the code is an averaged histogram 
combining 1000 outcomes of the simulation. For comparison, we generate the distribution 
of dg for the real AGN candidates in W5 and compare the results in Figure 14 (leftmost 
panel). In the right-hand panels we compare the output magnitude and color distributions 
of the simulated AGN (averaged over 1000 outcomes) with our W5 AGN candidates. 

These distributions can be used to generate a probability that an object in our AGN 
list, of given d^, [4.5] magnitude and [4.5] — [8.0] color is a YSO. In Table 6 we present the 
source properties of the 729 W5 AGN candidates. The full machine-readable table is given in 
the online material for this paper. Each object is listed with its coordinates and photometry 
in 2MASS and Spitzer hands, and given a combined 'P(YSO)' value, that is the probability 
given its local surface density, [4.5] magnitude and [4.5] — [8.0] color. 
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We statistically return a sample of AGN to the YSO list based on these probabihties to 
test the effect on the clustered fraction. We return samples based on the density probability 
alone, and then the combined density, color and magnitude probability. We run our MST 
clustering analysis on the resultant list of YSOs plus YSO- AGN. In Table 7 below we show 
the averaged results of adding AGN back to the YSO sample according to the different filters 
and compare to the original YSO-only sample result. 

The space density filter typically returns ^460 objects to the list. The combined filter 
returns ~230. Although a significant fraction of objects that are returned from AGN to YSOs 
join the clustered population, many are added to the distributed sources which increases the 
fraction outside groups and clusters. 

The objects returned to the stellar sample will also affect the Class I/Class II ratio. 
Applying the density filter returns a sample with ~16% Class I objects and the remainder 
Class II. If we apply the combined filter, typically 4% are Class I, 96% Class II. The effect 
on the ratio test is small: we increase the proportion of Class I objects in the Hll cavity 
from 4% to ~7%. The fraction associated with the molecular gas remains at 23%, thus our 
earlier result and conclusion still holds. A summary of the results of these results is given in 
Table 7 below. 

For the smaller sample of PAH galaxies, statistical analysis carried out as above pre- 
dicted on average 211 objects. However, we find 198 candidate PAH galaxies in W5, of which 
certainly some fraction are stellar objects as can be seen from the non-uniformity in Figure 
13. Although our model predicts an excess of PAH galaxies over that detected, we can still 
use the same procedure for these objects as for the AGN to generate a probability P(YSO), 
based on colors and space density. We present the 198 PAH galaxy candidates in Table 6 
along with the AGN with the same format of coordinates and photometry We provide a 
P(YSO) value for each, derived from local space density and [3.6] — [5.8] and [4.5] — [8.0] 
colors. 
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Note. — Table 6 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its 
form and content. Values in parentheses signify error in last 2 digits of magnitude value. Right ascension and Declination coordinates are J2000.0. 

''Galaxy type: AGN or 'PAH galaxy' as defined in text. 

''Probability that source is a YSO based on local space density, color, magnitude criteria: see text. 
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Fig. 13. — Distribution of extragalactic contaminants as classified by Spitzer photometry: 
AGN (blue points) and PAH galaxies (red). 



Table 7. Clusters/Groups in W5 





MST straight-line fit 
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% Clustered 
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66.0 0.8 


19 


43.4 


0.5 26 


Combined 


68.0 0.8 


17 


45.5 


0.5 25 


No AGN 


70.0 0.86 


16 


44.2 


0.54 27 
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Fig. 14. — Distributions of nearest neighbor distances, [4.5] magnitudes and [4.5] — [8.0] colors 
for objects classified as AGN in W5 (black histograms). Plotted in gray are the simulated 
distributions of simulated AGN according to our completeness estimates. We attribute the 
excess of objects in the W5 sample to young stars misclassified as AGN by our classification 
scheme. 
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